Charge Symmetry Violating Contributions to Neutrino Reactions 

J.T.LondergarQ 

Department of Physics and Nuclear Theory Center, 
Indiana University, 
Bloomington, IN 47405, USA 

A.W.ThomaJ] 

Department of Physics and Mathematical Physics, 
and Special Research Center for the Subatomic Structure of Matter, 
University of Adelaide, Adelaide 5005, Australia 
(Dated: February 7, 2008) 

Abstract 

The NuTeV group has measured charged and neutral current reactions for neutrinos on iron 
targets. Ratios of these cross sections provide an independent measurement of the Weinberg 
angle. The NuTeV value for sin 2 9 W is three standard deviations larger than the value measured 
in other electroweak processes. By reviewing theoretical estimates of parton charge symmetry 
violation (CSV), we study CSV contributions to the NuTeV measurement. We conclude that 
charge symmetry violating effects should remove roughly 30% of the discrepancy between the 
NuTeV result and other determinations of Sill U w • 
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I. INTRODUCTION 



In recent years precise experiments have provided detailed information regarding the 
parton structure of the nucleon. One of the more striking developments was the measurement 
of significant differences between up and down antiquark distributions in the nucleon sea. 
The first clear evidence for this was obtained from muon DIS on deuterium by the NMC 
group 111, l3l , that enabled a precise determination of the Gottfried Sum Rule 31. Later this 
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was more directly confirmed by Drell-Yan measurements in pp and pD reactions |4J, la, |6( 
and by semi- inclusive electroproduction at HERMES 7}. This sea quark flavor asymmetry, 
which had been anticipated [8j on the basis of spontaneously broken chiral symmetry 
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has been incorporated into the latest phenomenological parton distributions 

Another approximate symmetry in parton distributions is charge symmetry, which in- 
volves rotation of 180° about the two axis in isospin space, and which corresponds to inter- 
change of protons and neutrons (simultaneously interchanging up and down quarks). Our 
faith in the charge symmetry of parton distributions is justified from experience in nuclear 



physics, where this symmetry is respected to a high degree of precision |12|, [13J. Until re- 
cently, the quantitative evidence which could be extracted from high energy experiments, 
although not particularly precise, was consistent with charge symmetric parton distributions 
[l4 1 . As a result, all phenomenological analyses of deep inelastic scattering data in terms of 
parton distribution functions assume charge symmetry. Experimental verification of charge 
symmetry was difficult, first because high precision experiments were necessary to isolate 
charge symmetry violating (CSV) effects, second because the best tests required comparison 
of electromagnetic and neutrino-induced reactions, and third because CSV often mixes with 
quark flavor asymmetry effects. 

Recent precise experiments have now significantly decreased the upper limits on parton 
CSV contributions. The NMC measurements of muon DIS on deuterium [l| provide values 
for the cha rged lepton structure function FJfix, Q 2 ). In a similar Q 2 regime the CCFR Col- 
laboration [1^1 extracted the structure functions F^x^Q 2 ) from neutrino-induced charge- 
changing reactions. In sec. |H1 we review the comparison of these structure functions and 
we show that upper limits of a few percent can be placed on parton CSV in the region 
0.01 < x < 0.4. In sec. IIHI we review theoretical models for charge symmetry violation 
and we present two models for CSV parton distribution and an estimate of CSV parton 



distributions using phenomenological parton distributions. 

Recently, the NuTeV experimental group has measured both neutral and charged-current 
cross sections for neutrinos and antineutrinos on iron targets [16]. As was originally pointed 



out by Paschos and Wolfenstein 17(, ratios of these cross sections on isoscalar targets can 
provide an independent measurement of the Weinberg angle, sin 2 6^. The value of sin 2 ^ 
extracted by the NuTeV group is three standard deviations larger than the measured fit to 
other electroweak processes 18j . This striking result has led some people to suggest that 
they may be seeing evidence of physics beyond the Standard Model. 

Because of the importance of this result, it is crucial that all the effects contributing 
to this result be estimated to the best of our abilities. The NuTeV group has recently 
published a paper [19J that estimates the contributions to their result from three sources: 
the excess of neutrons in their iron target; the possible contribution from strange quarks and 
the effects of parton charge symmetry violation. As we will discuss in the following section, 
the isoscalar correction (excess neutrons in iron) makes a rather large, but apparently well- 
determined, correction to the extracted value of the Weinberg angle. At present, even the 
sign of the contribution from strange quarks is uncertain. In this paper, we will show that 
we can predict the sign of the contribution from parton CSV with some confidence. We will 
demonstrate that different theoretical predictions of the magnitude of the CSV correction 
to the neutrino determination of the Weinberg angle are in reasonable agreement with each 
other. CSV contributions should decrease the discrepancy between this experiment and the 
value of the Weinberg angle extracted from electroweak experiments in the vicinity of the 
Z mass. This will be reviewed in Sec. IIVI 

II. EXTRACTION OF WEINBERG ANGLE FROM NEUTRINO SCATTERING 

Paschos and Wolfenstein 17[ showed that one could obtain an independent measurement 
of the Weinberg angle, by taking ratios of charged-current and neutral-current cross sections 
for neutrinos and antineutrinos on isoscalar targets. They proposed measuring the ratio 

In Eq. (&hc) i s the neutral-current inclusive cross section, integrated over x and y, for 
neutrinos on an isoscalar tareret. The ouantitv On = M„J(M„ cos^l is one in the Standard 



Model. Alternatively, Eq. ^ can be written as 
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The NuTeV group used the Sign Selected Quadrupole Train beamline at FNAL to separate 

neutrinos and antineutrinos arising from pion and kaon decays following the interaction of 
800 GeV protons. The resulting interaction events were observed in the NuTeV detector, 
and were required to deposit between 20 GeV and 180 GeV in the calorimeter. CC and 
NC events were distinguished by the event length in the counters, as CC events contained 
a final muon that penetrated substantially farther than the hadron shower. 

Rather than working directly with the Paschos-Wolfenstein ratio, the NuTeV collabora- 
tion measured the individual ratios R v and R v defined in Eq. |2J and took the value of r 
from earlier experiments. From the result, R v = 0.3916 ± 0.0007 and R v = 0.4050 ± 0.0016, 
they extracted sin 2 6 W = 0.2277 ±0.0013 (stat) ±0.0009 (syst). This value is three standard 
deviations above the measured fit to other electroweak processes, sin 2 ^ = 0.2227± 0.00037 
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|18| . This result is striking, and if no other effects can account for this discrepancy, it may 
be evidence of physics beyond the Standard Model. 

Several approximations have been made in deriving Eq. It is true only for isoscalar 
targets, includes only the contributions from light quarks, and assumes the validity of parton 
charge symmetry. The NuTeV group has recently investigated how their result changes when 
these assumptions are removed |l9j . The corrections to the Paschos-Wolfenstein ratio of Eq. 
Qtake the form 



AR- 



U v + D, 



2A 2 + 3(A" + A 



(3A 2 + A 2 ) 
(U v + D Y ) 

Q = xq(x)dx 
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-5N (U v ~D V ) + ^ (5U V - 5D V 



Q v = Q-Q (3) 

In Eq. El the quantity Q is the total momentum carried by a quark of flavor q, and the 
quantity Q v is the total momentum carried by valence quarks of that flavor. 5N = (N—Z)/A 
is the fractional neutron excess, A 2 d = (e^ d ) 2 — (e^) 2 and e c is the kinematic suppression 
factor for massive charm nroduction. 



The isoscalar contribution (the term proportional to SN in Eq. EJ) is straightforward to 
calculate. Using the values 5N = .0567, 3A 2 + A 2 , = .4804, and {U v - D v ) / (U v + D v ) ~ 0.46, 
gives an isoscalar correction to sin 2 6 W of about —0.0125. The NuTeV group reports an 
isoscalar correction of —0.0080, with a very small error [20]. This differs from the 'naive' 
correction of Eq. |H] because the NuTeV group corrects for factors like experimental cuts, 
experimental backgrounds and the enhanced sensitivity of their experiment to neutrino scat- 
tering at low x. All of these factors were used as input in a detailed Monte Carlo simulation 



of their experiment. However, Kulagin |2l( has recently claimed that the uncertainties in 
the isoscalar corrections are likely to be considerably larger than estimated by the NuTeV 
group. 

The contribution from strange quarks depends on the quantity S v = S — S, the difference 
between the momentum carried by strange quarks and strange antiquarks. Even the sign 
of this quantity is not firmly established. Barone et al. [22] analyzed the CDHS neutrino 
charged- current inclusive cross sections and charged lepton structure functions j^J, and 
argued that some improvement is obtained by_ allowing an asymmetric strange sea with 
S v > 0. However, the CCFR [^J and NuTeV pij charged-current and dimuon results show 
significant disagreement with the CDHS results at large x. The NuTeV group finds a best 
fit to their results with a slightly negative value S v = —.0027 ± 0.0013. 

In the remaining sections, we will review the origin of valence quark charge symmetry 
violation, and we will apply several theoretical models to estimate the CSV contribution to 
the NuTeV value for sin 2 9 W . 

III. VALENCE QUARK CHARGE SYMMETRY VIOLATION 

Because the Paschos-Wolfenstein relation involves the difference between neutrino and 
antineutrino cross sections, it is dominated by contributions from valence parton distribu- 
tions, 

q v (x) = q(x) -q(x) , (4) 

for a particular quark flavor q. We can gain insight into the origin and magnitude of parton 
charge symmetry violation by using a method for calculating twist-two valence parton distri- 
butions developed by the Adelaide group 3 . This method evaluates quark distributions 



through the relation 

q(x^ 2 ) = M £ I W + (0)|iV>| 2 5(M(l " x)-p$) (5) 
x 

In Eq. ip + = (1 + a^)ip/2 is the light front operator that removes a quark or adds an 
antiquark to the nucleon state \N), /i 2 represents the starting scale for the quark distribution, 
\X) are all possible final states that can be reached with this operator, and p x is the plus 
component (p + = p 3 + E(p)) of the residual system. Therefore, |X) = 2q, 3q + q,4q + 2q, . . .. 

Thomas and collaborators showed that one could obtain quark distributions in semi- 
quantitative agreement with experiment using simple quark models such as the MIT bag 
27| in Eq. then evolving the resulting quark distribution from the bag scale fi 2 to the final 
value of Q 2 . The contribution to a quark distribution from intermediate state X produces 
a peak at a value x p ~ 1 — M x /M, where M x is the effective mass of the state X. We 
note that for two-quark intermediate states, x p ~ 1/3, while x p is negative for states with 
four or more quarks. Consequently, for x > 0.2, where valence quarks dominate, qualitative 
estimates can be obtained by including only two-quark intermediate states in Eq. Most 
importantly, this observation means that the major qualitative features of the spin and flavor 



dependence of the valence parton distribution functions can be understood simp 
of the hyperfine mass splitting between two-quark states with spin zero or one 
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similar analysis can also explain the spin and flavor dependence of strange baryon parton 



distributions at large , M , and also q uark fragmentation tactions at large energy fraction 
z |30j. 

We can use Eq. 03 to estimate charge symmetry violating effects, e.g., the difference be- 
tween the up quark valence distribution in the proton and the down quark in the neutron, 
Sd v (x) = d,P(x) — Uy(x). There are four sources of CSV contributions: charge symmetry 
violation in the quark wavefunctions; electromagnetic effects that break charge symmetry; 
mass differences of the struck quark; and mass differences in the spectator multiquark sys- 
tems. Model quark wavefunctions are found to be almost invariant under the small mass 
changes typical of CSV. At sufficiently high energies, electromagnetic effects should also be 
small, and these are neglected. Consequently, parton charge symmetry violation will arise 
predominantly through mass differences rrid — m u of the struck quark, and from mass differ- 
ences in the spectator multi-quark system. Both of these contributions will result in small 
shifts in the argument of the enerffv-conservinGr delta, function in Eo.lKl Since at laxee x the 



contribution to the parton distribution is dominated by the two-quark intermediate state, 
we can thus obtain quantitative estimates of the sign and magnitude of parton CSV from 
these terms. We stress that the change in the mass of the spectator pair is exactly the same 
mechanism that leads (through a much larger mass difference) to an understanding of the 
major features of the spin and flavor dependence of valence distributions, therefore one has 
a fair degree of confidence in this particular term. 
Consider the "minority valence quark" CSV term, 

6d v (x) = d p Y (x) - <(x) . (6) 

The main contribution to minority valence quark CSV arises from the mass of the diquark 
system, which is uu for the proton and dd for the neutron. The two-quark contribution to 
the valence quark distribution will peak at x v v ~ 1 — M uu /M for d v and x™ ~ 1 — Mdd/M 
for Uy, thus the down quark distribution in the proton will be shifted to higher x and the 
up quark distribution in the neutron will be shifted to lower x. This means that, at large x, 
5d v (x) will be positive. 

Next, we consider the "majority valence quark" CSV term, 

5u v {x) = - <0) . (7) 

The mass of the diquark state is ud for both and d™. As a result, the two-quark contri- 
bution to the majority quark valence distribution will peak at x v v ~ 1 - M ud /MP for up and 
Xp ~ 1 — M u d/M n for d™. As a result we expect 5u v (x) to be negative at large x. 

An important constraint on CSV parton distributions is that they respect the normal- 
ization of valence quarks. The integral over x of valence quark distributions must give the 
total number of valence quarks, i.e., 

u p Jx)dx = f 1 d™{x)dx = 2 
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d p v (x)dx = [ u™(x)dx = 1 , (8) 
o Jo 

and hence Eq. |H] requires that 

/ 5d v (x)dx= / 5u v (x) dx = 0. (9) 
Jo Jo 

It is important that valence quark CSV distributions obey this constraint, otherwise one is 
effectively chanffinGr the total number of ud or down valence nuarks in the nucleon. 
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FIG. 1: Valence quark CSV contributions, xSq v (x) vs. x. Solid line: x5u v ; dash-dot line: x5d v . 
Calculated using MIT bag model wavefunctions by Rodionov et al. , Ref. |3l(, and evolved to 
Q 2 = 10 GeV 2 . 

Since we can infer the sign and relative magnitude of the CSV parton distributions at 
large x, and these distributions must integrate to zero, we can obtain at least qualitative 
values for the valence CSV distributions. 6d v (x) will be positive at large x and thus negative 
at small x, while Su v (x) will be negative at large x and positive at small x. In Fig. Q we 
show x5q v (x) for valence up and down distributions as calculated by Rodionov, Thomas 
and Londergan [3.1]. These distributions are calculated from Eq. using an MIT bag model 
for the quark wavefunctions. These numerical calculations included all four sources of CSV 
noted earlier. The resulting CSV distributions were evaluated at the bag scale and evolved 
upwards in Q 2 . The distributions shown in Fig. ^ are evolved to Q 2 = 10 GeV 2 . The parton 
distributions exhibit the qualitative effects that we inferred; at large x, 5d v (x) > and 
Su v (x) < 0. The CSV distributions must change sign at small x to satisfy the requirement 
of valence quark conservation, summarized in Eq. EH 

Charge symmetry violating parton distributions were also calculated by Sather j^. 
Starting also from Eq. Sather derived the following analytic approximation relating CSV 




distributions to valence parton distributions, 

5d v ( x ) = -^KWl-^iW 

5u Y (x) = *Jt (-± [xu v (x)} + ^ v (*)j (10) 

In Eq. 8M = M n — M p = 1.3 MeV is the neutron-proton mass difference, 8m = rrid — m u 
is the down-up quark mass difference, and M is the average nucleon mass. To calculate 
the CSV parton distributions, Sather used parton distributions obtained from the MIT bag 
model in Eq. [TU| Qualitatively, Sather's CSV parton distributions are quite similar to those 
of Rodionov et al. jsi|. However, 8u v (x) is smaller than the corresponding quantity for 
Rodionov et al. , and Sather's CSV distributions peak at somewhat smaller x. Sather's 
CSV distributions are qualitatively similar to those shown in Fig. E3 

Both of these charge symmetry violating distributions were constructed using bag model 
wavef unctions. Such wavefunctions can give qualitative agreement with phenomenological 
parton distributions, but these model-generated parton distributions systematically under- 
predict the results of phenomenological distributions at large x. Here, we present a method 
for generating valence quark CSV parton distributions from phenomenological distributions. 
We will use the approximate formula derived by Sather j^J , Eq. El that relates CSV parton 
distributions to the derivatives of valence quark parton distributions. 

Although Sather obtained his CSV distributions by using parton distributions from the 
MIT bag model in Eq. El we will insert phenomenological parton distributions into this 
equation. There is a problem with this approach. As we pointed out in Eq. conservation 
of valence quark probability requires that the integral over all x of Sq v in Eq. El be zero. 
Since the terms in this equation are just derivatives of parton distributions, integration over 
x simply involves evaluating the parton distributions at zero and one. However, phenomeno- 
logical parton distributions go like x~ 1 / 2 in the limit x — > 0, hence when phenomenological 
distributions are used the integral of the CSV distributions in Eq. El will not be zero, in 
fact the integrals will blow up. 

The problem originates because Eq. El is a reasonable approximation for the parton 
distribution that arises when a nucleon consisting of three valence quarks splits into a quark 
and diquark. This gives the dominant contribution to the parton distribution at large x. 
However, at small x the valence distribution is dominated by higher mass Fock states that 
include manv ouark-antinuark nairs. For states involving such larere excitations, the effects 



TABLE I: Parameters for valence quark distribution from CTEQ4LQ parton distribution, Ref. 
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1.315 


-0.427 


3.281 


10.614 


0.607 




0.852 


-0.427 


4.060 


4.852 


0.266 



of quark and nucleon mass differences should be negligible. However, Eq. El incorrectly 
predicts very large contributions at small x. Consequently, we need to suppress the large 
CSV effects produced by Eq. ED at small x. We will deal with this problem in a simple, 
and completely phenomenological, way. Parton distributions from the CTEQ collaboration 
Q 



351 ] have been parameterized using the form 

q v (x)=N[x a + Cx~<}(l-x) p (11) 



The small- a; behavior is governed by the parameter a ~ —0.5. We want a modified parton 
distribution that will vanish at small x, leaving the large-x behavior relatively unchanged. 
We thus replace the CTEQ parton distributions q v (x) in Eq. HHwith q v (x), defined by 

q v (x) = N[x a + Cx^\ [(1 - xf - (1 - xf +12 ] (12) 

By inspection, q v (x) vanishes at x — 0, and for large x there is a very small difference between 
the modified parton distribution of Eq. El and the phenomenological parton distribution. 

We calculated valence quark CSV distributions using Eq. HI! with the CTEQ4LQ phe- 
nomenological parton distributions j^], modified using Eq. ^] The coefficients for the 
CTEQ4LQ parton distributions, appropriate for Q 2 = 0.49 GeV 2 , are given in Table [I] 

Sather's analytic expression, Eq. [TH1 is appropriate for a nucleon at a low starting scale, 
of order Q 2 ~ 0.25 — 0.5 GeV 2 . The CTEQ4LQ parton distribution was introduced 
specifically for starting scale Q 2 = 0.49 GeV 2 . The resulting CSV distributions were evolved 
to the higher value Q 2 = 20 GeV 2 appropriate for the NuTeV experiment, using the QCD 



evolution program of Miyama and Kumano 



The evolved CSV parton distributions are 



shown in Fig. |2J The solid curve shows x5d v and the dash-dotted curve gives xSu v . 
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FIG. 2: Valence quark CSV contributions, x5q v (x) vs. x. Solid line: x5d v ; dash-dotted line: 
x5u v . Calculated from valence quark distributions using the analytic approximation of Sather, 
Ref. U, and the CTEQ4LQ parton distributions of Ref. [3J], modified according to Eq. 1121 
These distributions were then evolved from the starting scale Q 2 = 0.49 GeV 2 to Q 2 = 20 GeV 2 
appropriate for the NuTeV experiment. 

IV. CHARGE SYMMETRY VIOLATING CONTRIBUTIONS TO NEUTRINO 
REACTIONS 

Using the valence quark CSV distributions reviewed in Sec. Hill we can estimate the CSV 
contribution to the extracted value of sin 2 ^. From Eq. El the CSV corrections to the 
Paschos-Wolfenstein ratio are of the form 

Q v — x q v (x) dx 
Jo 

5D V = x \d p v (x) - u^ix)] dx 
Jo 

5U V = [ x [vPJx] - <%(x)) dx (13) 
Jo 

In Table |H] we show the contributions of the different theoretical CSV estimates to the 
Paschos-Wolfenstein ratio. We have broken down the individual contributions so that AU V 
and AD V are the total CSV effects arising from SU V and SD V , respectively. From Fig. [I\ 
121 in all cases both 5U V and 5D V make a negative contribution. Therefore the net CSV 
contribution will be negative. This will decrease the discrepancy between the value of 
sin 2 6 W extracted in the NuTeV experiment, and the best value obtained from high-energy 
electroweak interactions. 



TABLE II: CSV corrections to determination of sin 2 9 W in neutrino scattering. PW is contribution 
to the Paschos-Wolfenstein ratio, Nu is weighted by the NuTeV functional. AU (AD) is total 
contribution from 5u v (5d v ), and Tot is total CSV correction. 
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AD Nu 
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Rodionov 


-.0010 


-.0011 


-.0020 


-.00065 


-.00081 


-.0015 


Sather 


-.00078 


-.0013 


-.0021 


-.00060 


-.0011 


-.0017 


analytic 


-.00075 


-.0013 


-.0021 


-.0005 


-.0009 


-.0014 



Table |H] shows that CSV contributions to the Paschos-Wolfenstein ratio from three the- 
oretical models give almost identical results, —0.0020 or —.0021. Since the value of sin 2 ^ 
from the NuTeV experiment, without CSV corrections, was 0.005 larger than the best value 
from electromagnetic interactions, our calculated CSV effect would reduce the discrepancy 
between the neutrino and electromagnetic measurements of sin 2 9 W by 40%. 

However, as pointed out by the NuTeV group [19], it is not appropriate to use Eq. [LUto 
determine the CSV effects on the value of sin 2 9 W , as the NuTeV extraction of this quantity 
does not rely on the Paschos-Wolfenstein ratio, but instead uses the absolute ratios R u and 
R v defined in Eq. |2J and compares the results with a full Monte Carlo simulation of the 
experimental processes. The NuTeV group produced functionals that give the sensitivity of 
their observables to various effects. These are summarized in a single integral, 

AS = f 1 F[S,5;x] 5{x)dx (14) 
Jo 

Eq. El gives the change in the extracted quantity £ resulting from the symmetry violating 
quantity 5(x). The functionals appropriate for the observable sin 2 9 W and the parton CSV 
distributions, were provided in Ref. 19j. 

We have taken the parton CSV distributions and folded them with the NuTeV functionals. 
The net CSV correction to the value of sin 2 6 W is listed as Tot^u in Table ILT1 The CSV 
contributions are still negative, i.e., they decrease the discrepancy between the neutrino 
and electromagnetic measurements of the Weinberg angle. However, the CSV change in 
sin 2 6 W is slightly smaller than estimated in Eq. ED This is because the NuTeV experiment 
is somewhat more sensitive to small-x physics than to effects at larger x. The CSV parton 
distributions, weighted by x, change sign at small x and reach a maximum at larger x. 



Thus the functional, which emphasize the CSV distributions at small x where they are 
small and change sign, and de-emphasize them at larger x, give a smaller CSV effect than 
predicted by the Paschos-Wolfenstein ratio. Nevertheless, the CSV contributions to the 
NuTeV result, weighted with their functionals, range from —0.0014 to —.0017, and thus 
reduce the discrepancy between the neutrino and electromagnetic measurements of sin 2 9 W 
by roughly 30%. 

The NuTeV group estimated the CSV correction and obtained a much smaller value 
than ours. However, in order to obtain CSV parton distributions, they took the ratio 
5q v (x)/q v (x) from Rodionov et al. jjilj], and multiplied this ratio by parton distributions 
determined from neutrino scattering. Since the ratio was determined using parton distribu- 
tions from a quark model, and the parton distributions were obtained from quite a different 
source, we are not convinced of the accuracy of the resulting CSV distributions. In partic- 
ular, CSV distributions constructed in this way will not satisfy any relation such as Eq. ITUl 
nor will they satisfy valence quark conservation of Eq. H3 



V. CONCLUSIONS 



We have reviewed the corrections that parton charge symmetry violation should make 
in determining sin 2 9 W in neutrino scattering. Although parton CSV effects are sufficiently 
small that neither their sign nor magnitude has been measured to date, we argued that at 
large x the CSV distribution Sd v (x) should be positive and 5u Y (x) should be negative. The 
CSV distributions must preserve the overall number of valence quarks in the neutron, hence 
from Eq. El both CSV distributions must change sign at small x. 

We investigated two theoretical models of parton CSV. In both cases, the parton distribu- 
tions are calculated from a model of QCD, such as the MIT bag, and the models determined 
how the CSV distributions could be related to the parton distributions and their deriva- 
tives. In a third case, we inserted phenomenological parton distributions from the CTEQ 
group into analytic expressions relating the CSV distributions to parton distributions. We 
removed the tendency of these analytic expressions to (incorrectly) give too much weight to 
the small-x region, by damping the phenomenological distributions at very small x. In all 
cases, the CSV distribution was obtained at a low starting scale, and then evolved to higher 
Q 2 , more appropriate for the Q 2 value of the neutrino experiments. 



All of these theoretical CSV distributions serve to decrease the value of sin 2 6 W extracted 
from the NuTeV experiment; the size of these corrections was remarkably similar, ranging 
from —.0014 to —.0017. Since the value of sin 2 9 W extracted from neutrino scattering was 
greater than that obtained from electromagnetic interactions by +.005, inclusion of a charge 
symmetry violating contribution of this magnitude would reduce this discrepancy by about 
30%. We emphasize that CSV effects have yet to be confirmed by direct experiment, so our 
calculations only give the best theoretical estimate of these contributions. The predicted 
CSV effects arise from the same mechanism that correctly predicts the spin and flavor 
dependence of valence quark distributions, so we would be surprised if our CSV effects did 
not have the correct sign and roughly the right magnitude. It is rather remarkable that 
although the predicted corrections from parton charge symmetry violation are quite small, 
high energy experiments have now reached a precision where even small absolute effects have 
a significant impact on our ability to extract fundamental quantities. 
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